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^ Abstract 

I Is it possible to infer the time evolving quantum state of a multichromoplioric system from a sequence 

of two-dimensional electronic spectra (2D-ES) as a function of waiting time? Here we provide a positive 

answer for a tractable model system: a coupled dimer. After exhaustively enumerating the Liouville 

(— I , pathways associated to each peak in the 2D-ES, we argue that by judiciously combining the information 

^ \ from a series of experiments varying the polarization and frequency components of the pulses, detailed 
a ■ 

■ information at the amplitude level about the input and output quantum states at the waiting time can 
O^. be obtained. This possibility yields a quantum process tomography (QPT) of the single-exciton manifold, 
which completely characterizes the open quantum system dynamics through the reconstruction of the 

> 

\Q ' process matrix. This is the first of a series of two articles. In this manuscript, we specialize our results 
' to the case of a liomodimer, where we prove that signals stemming from coherence to population transfer 
and viceversa vanish upon isotropic averaging, and therefore, only a partial QPT is possible in this case. 
However, this fact simplifies the spectra, and it follows that only two polarization controlled experiments 
(and no pulse-shaping requirements) suffice to yield the elements of the process matrix which survive under 
^ . isotropic averaging. The angle between the two site transition dipole moments is self-consistently obtained 

H ; 

, from the 2D-ES. Model calculations are presented, as well as an error analysis in terms of the angle between 
the dipoles and peak overlap. In the second article accompanying this study, we numerically exemplify 
the theory for heterodimers and carry out a detailed error analysis for such case. This investigation 
provides an important benchmark for more complex proposals of quantum process tomography (QPT) via 
multidimensional spectroscopic experiments. 
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Multidimensional optical spectroscopies (MDOS) provide very powerful tools to study excited 
state dynamics of multichromophoric systems in condensed phases. These techniques distribute 
spectral features along several dimensions, uncluttering data which would otherwise appear ob- 
scured in linear spectroscopies and simultaneously yielding novel information on the dynamics of 
the probed system Ipj]. Possibilities in multidimensional techniques include decongesting spectral 
lineshapes, differentiating between homogeneous and inhomogeneous broadening mechanisms, pro- 
viding unambiguous signatures about couplings between chromophores, and yielding signatures of 
coherent and incoherent processes involving excited states at the amplitude level 2|, |3|. Although 
MDOS have been historically inspired by their NMR analogues, the timescales of the physical and 
chemical processes studied through MDOS are quite different from the ones in NMR [#i7|. The 
characteristic timescales of NMR are milliseconds, a resolution that does not allow for the obser- 
vation of a wide variety of chemical dynamics in condensed phases ocurring in the orders of femto 
and picoseconds. On the other hand, femtosecond timescales can be easily accessed with ultrafast 
optical techniques. Examples of phenomena studied via MDOS are vast and include molecular 



reorientation processes and solvation dynamics jS 
in organometallic complexes 



systems 



18 



11 



9|, electron transfer 



13| or halogens in rare gas matrices 



in carbon nanotubes [16], protein unfoldin g ki netics 



quantum wells 



14, 



vibrational coherences 



15] . phonon dynamics 



17] . excitonic dynamics in light-harvesting 



23] and organic polymers 24 J25], as well as many-body physics in semiconductor 



261-128] and quantum dots 
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Traditionally, the spectroscopy of condensed phases is formulated as a response problem: The 
molecular system is perturbed with a sequence of short laser pulses and the coherent polarization 
response due to this set of perturbations (nonlinear polarization) is subsequently measured [l]. In- 
formally, we can describe the exercise as 'kicking' the quantum black box (molecular system) and 
'listening to the whispers' (measuring the response) due to the kicks, from which some properties of 
the box can be inferred. This description of spectroscopy is reminiscent to an idea stemming from 
the quantum optics and quantum information processing (QIP) communities, namely, quantum 
process tomography (QPT) 30M33]. Broadly speaking, QPT is a systematic procedure to char- 
acterize a quantum black box by sending a set of inputs, measuring their outputs, and analyzing 
the functional relationships between them. With the increasing effort of quantum engineering of 
gates and devices, QPT constitutes a cornerstone of QIP theory and experiment, as it provides a 
necessary check on the performance of the respective quantum black boxes. A natural question 
arises from the comparison of the two aforementioned concepts: Can the spectroscopy of condensed 
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phases be formulated as a QPT? In a previous study j34| . we provided an affirmative answer to this 
question, at least for a molecular dimer. We showed that a set of two-color polarization controlled 
rephasing photon echo experiments is sufficient to reconstruct the density matrix elements associ- 
ated with the dynamics of the single exciton manifold, and therefore, systematically characterize 
the excited state dynamics of the dimer, which can be regarded as the black box. For pedagogical 
reasons, we found it simpler and more convenient to concentrate our attention in the real time 
picture of the experiment, to make an explicit identification of the preparation, evolution, and 
detection steps of the QPT, with the coherence, waiting, and echo times, respectively. However, 
due to the widespread practice of displaying partially Fourier transformed data of the nonlinear 
optical polarization with respect to certain time intervals, it is worthwhile translating our results 
to the more visual two-dimensional electronic spectrum (2D-ES), and in fact, this is one of the 
main results of this work. 

The present article is organized as follows: We begin in section I with a review of some relevant 
ideas of QPT and also introduce the process matrix as the main object to be reconstructed by 
means of QPT. In section II, relevant details on the dimer model system are presented. Section 
III describes the rephasing heterodyne photon echo experiment for the dimer and explains that 
the collected macroscopic polarization signal is a linear combination of elements of the process 
matrix at the waiting time x{T)- This implies that QPT can be performed by repeating several 
experiments with different pulse parameters in order to extract these elements. In section IV, the 
ideas of section III are mapped into the language of a 2D-ES, where each of the diagonal and 
cross peaks is associated with a set of elements of x{T), and each of the axes of the spectrum 
can be associated with a preparation and a detection stage. Finally, in section V, these ideas are 
specialized to homodimer systems where, after isotropic averaging, only a partial QPT is possible, 
as some elements of x{T) are undetectable. Nevertheless, we note that the partial QPT is easily 
realized with current experimental capabilities, since it can be reconstructed with only two spectra 
resulting from different pulse polarization configurations for each given waiting time. The angle 
between the site dipoles is self-consistently obtained from these spectra, and an error analysis based 
on this angle as well as peak overlaps is carried out. Numerical calculations on a secular Redfield 
dissipation model are presented. A detailed analysis for heterodimers is carried out in the next 
article accompanying this investigation. Extensions of the procedure to account for inhomogeneous 
broadening, more sophisticated signal analysis, as well as bigger systems, are discuss at the end of 
this manuscript. 
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I. RELEVANT CONCEPTS OF QUANTUM PROCESS TOMOGRAPHY AND GEN- 
ERAL DEFINITIONS 



Consider an arbitrary quantum system (quantum black box) interacting with an environment. 
We are interested in its evolution as a function of time T in the form of a reduced density matrix 
p(T). Very generally, this evolution is a linear transformation on the initial quantum state jsS, 36 |: 



p(T)=x(T)p(0). 



(1) 



x{T) is the central object of this article, and shal 
as an integrated equation of motion for every T 
for the Liouville space of the system: 



be called process matrix. Eq. ([T]) can be regarded 



10l|. Eq. ([T]) can be expressed in terms of a basis 



Pab{T) = J2Xabcd{T)p,d{0). (2) 

cd 

For purposes of this article, we present two useful definitions. Consider the Liouville space C 
of the system, and classify the vectors of C in proper and improper density matrices. A state or a 
density matrix is proper if it satisfies all the conditions of a physical quantum state; namely, this is 
Hermitian, positive semidefinite, and has trace one. An improper state is any matrix that lives in 
the Liouville space but is not a proper density matrix. Clearly, any improper density matrix in the 
same Liouville space may be written as a unique linear combination of proper density matrices. 
In principle, Eq. ([2]), being a physical equation of motion, is restricted to the domain of proper 
density matrices p(0). However, by linearity, its extension to any linear combination of proper 
states is well defined, so its validity for improper density matrices is not under question. 

The meaning of the process matrix x(T) is easy to grasp: Conditional on the initially state 
being prepared at p(0) = |c)((i|, Xabcd{T) is the value of the entry ab of the quantum state after 
time T, p(T), i.e. Xabcd{T) = {a\p{T)\b). Therefore, Xabcd{T) denotes a state to state transfer 
amplitude. Note that p(0) = \c){d\ is an improper density matrix if c ^ d (coherences on their 
own are not valid quantum states). However, improper states are not necessarily unphysical as 
one expects at a first glance. Most of our intuition for nonlinear spectroscopies in the perturbative 
regime stems from the consideration of how a perturbative amplitude created at a certain entry 
\c){d\ of the total (proper) density matrix is transferred to other entries due to free evolution, as 
time progresses [1]. It is not the evolution of the total density matrix (which to leading order 
is unperturbed, mostly in its ground state, and not yielding a time-dependent dipole) what is 
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effectively monitored in tlie pliase-matclied signal, but the evolution of an effective density matrix, 
such as |c)((i|, which can be improper. Terms such as transfer from population to population, 
coherence to coherence, population to coherence, and coherence to population are all ubiquitous in 
the jargon of MDOS. However, the monitoring of the latter is often ambiguous, incomplete, and in 
most cases, qualitative. Obtaining quantitative information about these events amounts to finding 
each of the elements of x(T). 

The transformation in Eq. ([2]) is limited by two classes of restrictions for the process matrix 
associated with Hermiticity and trace preservation: 



Jed, 



^Xaacd{T) = 6c 

a 

Xabcd{T) = XbadciT). 



(3) 
(4) 



We derive these conditions in Appendix A, but their content is intuitive: If p(0) is a proper density 
matrix, p{t) remains as a valid quantum state as T evolves if these two requirements are preserved. 
In particular, elements of the form Xaabb{T), which denote population transfers, are purely real as 
one expects, whereas the other elements are in general complex. As a comment to our previous 
discussion, by linearity, these conditions must also be satisfied even if p(0) is improper (notice that 
x{T) does not depend on p(0)). 

Equations ([1]) and ([2]) are remarkable because they guarantee that, in principle, if x{T) is known, 
the quantum black box described by p(T) is perfectly understood, as it predicts by linearity the 
evolution of an arbitrary initial state of C. Although p{t) describes an open quantum system, 
details about the environment evolution need not be included explicitly, but only in an averaged 
sense in the elements of x{T)- We shall operationally define QPT as any procedure to reconstruct 
x{T)- A possible QPT is the following: (a) Prepare a linearly independent set of states p(0) 
that spans C; (b) for each of the prepared states, wait for a free evolution time T and determine 
the density matrix at that time. Any protocol for determining a density matrix for a system is 
called Quantum State Tomography (QST) 37H4l|. In essence, QPT can be carried out for any 
system if both a selective preparation of initial states and QST can be achieved. Variants of this 



methodology exist although all of them operate within the same spirit 



QPT has been 



successfully implemented in a wide variety of experimental scenarios, including nuclear magnetic 



resonance 



42|-|44| . ion traps j45[. single photons 46|, 4^, solid state qubits 48||, optical lattices 



49|, and Josephson junctions 



50] . In this article, we show how to perform QPT for a model 
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coupled heterodimer using two-color polarization controlled heterodyne photon-echo experiments, 
extending the domain of application of QPT to systems of chemical and biophysical interest. 

II. MODEL SYSTEM: COUPLED DIMER 



Consider a molecular dimer described by the effective Hamiltonian 




Hs = ujAa\aA + Wsa^a^ + J{a\aB + a%aA), (5) 

where af and are creation and anhilation operators for a single Frenkel exciton in the site 
i G {A,B}, ua,ujb are the first and second site energies, and J is the coupling between the 
chromophores. 

The standard diagonalization of this Hamiltonian, which is effectively a two-level system for the 
single-exciton manifold, follows from defining some convenient parameters: The average of the site 
energies co = ^{uja + (^b), the difference A = — (^s), and the mixing angle = ^ arctan (■^). 
By introducing the operators: 



tta = cos 6a A + sin 6a b , 

a/s = — sin 6aA + cos 6aB, (6) 
the Hamiltonian in Equation ([5]) can be readily written as: 

Hs = Uaa^aa + W/ja^a/j, (7) 
where the eigenvalues Ua and a;/3 of the single excitons are: 



cJq, = w + A sec 26, 

Ufj = u - Asec26. (8) 

Denoting \g) as the molecular ground state or the excitonic vaccuum, 1^4) = a^g) and \B) = a'^\g) 
are the excitons at each site, whereas \a) = a^\g), = a~^\g) are the delocalized excitons. The 
biexcitonic state, expressed by |/) = a\a~^\g) = c'^c^\g), also plays a role in our study, as it is 
resonantly accesed through excited state absorption (ESA) after several pulses. Notice that the 
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Hamiltonian Hs does not contain two-body operators, and therefore does not account for exciton- 

exciton binding or repulsion terms, so the energy level of the biexciton is just the sum of the two 

n 

exciton energies, w/ = Wq, + = + [2]. Defining the following relations hold: 



UJ^g = Uf^. (9) 

Since we are concerned with the interaction of the chromophores with electromagnetic radiation, 
we make some remarks on the geometry of the transition dipoles (see Fig. 1). Let = (z|/it|j). 
Assume that the transition dipole moments from the ground to the single excitons in the site basis 
are Hgj^ = fij^g = and /x^^ = ^^g = ds, respectively. It follows that the dipole moments 
for i,j G {a,/3,f} are located in the same plane, but in general have different magnitudes and 
directions: 



t^ag 




COS 6 sin 6 




dA 






— sin 9 cos 9 




ds 






sin 9 cos 9 




dA 






cos 9 — sin 9 




ds 



(10) 



In this model, we shall consider /Xj^ = /x^j. As enumerated in our model, dipole mediated transitions 
only couple the ground state to the single excitons, and the single excitons to the biexciton. 



III. PHOTON-ECHO EXPERIMENT AS QUANTUM PROCESS TOMOGRAPHY 

Consider a four-wave mixing experiment where an ensemble of identical dimers interacts with 
a series of three ultrashort laser pulses. The perturbation due to these pulses is given by: 

3 

V{r, t) = -A ^ /X ■ eiE{t - t.)e^*'«— + c.c, (11) 

i=l 

where A is the intensity of the elect ric fie ld, which is assumed to be weak, /x is the dipole operator. 



Cj, tj, fcj, uji denote the polarization 102|, time center, wavector, and carrier frequency of the i — th 
pulse, and r is the position of the center of mass of the molecule. E{t) is the slowly varying in 
time pulse envelope, which we choose as a Gaussian with width a, or full-width half-maximum 
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FWHM = 2 a/2 In 2a, E{t) = e-*'/(2'^'). The pulses are sent to the sample in a non-collinear fashion 
to the sample, generating a time-dependent dipole in each of the molecules. Since the characteristic 
size of a molecule is much smaller than the wavelength of the radiation, 27r/|fci|, each molecule 
only experiences a potential that changes in time but is uniform in space, in consistency with the 
dipole approximation. Nonetheless, the spatial dependence of the pulses is still important, as the 
phases e^**'" '^ are imprinted to molecules located across different positions r in the sample. The 
size of the sample is much larger than 2n/\ki\, so there is a considerable spacial modulation of the 
polarization due to these phases. Denoting the time-dependent state of the molecule at position 
r by p{r,t), a perturbative treatment allows us to decompose the density matrix into Fourier 
components: 



p(r,t) = 5^p,(t)e^*'-^ (12) 

s 

where kg = Iki + mfc2 + nk^ and /, n, m are integer numbers. Notice that kg equals to a linear 
combination of the wavevectors associated with each pulse. As expected, the action of a pulse 
on each molecule attaches a spatial phase to its quantum state, so the total phase accumulated 
by it equals e*'"*'^ for each combination of perturbations. Each improper density matrix Ps{t) 
corresponds to one of these phases, and can be calculated by keeping track of the actions of the 
pulses in the bra and the ket of the system using double-sided Feynman diagrams. Eq. (fT2|) 
implies that the optical polarization induced on the molecule can also be Fourier decomposed into 
different components 53|-|55|: P{r,t) = Tr{jl{r)p{r,t)) = Ps(t)e**'*"^, where fi{r) denotes the 
dipole operator of the molecule located at r. The experimental setting we describe is analogous 
to the one of an array of dipole antennas which are spatially phased in a grating with respect to 



56|. 



each other and oscillate in time. Classical electromagnetism predicts that the induced macroscopic 
polarization of this array emits radiation which is precisely concentrated along the vectors kg. This 
condition, which reflects conservation of momentum of the fields, is known as phase-matching 
A fourth pulse of the same wavevector as one of the kg, known as the local oscillator, is allowed 
to interfere with the radiation along that direction. By varying the phases of this fourth field, 
two heterodyne detections can be carried out to extract the real and imaginary components of 
Ps{t) = Ps{t) ■ 64 respectively, where 64 is the polarization of the local oscillator 

In this article, we are interested in the signal along kpE = —ki + ^2 + ^3, the so called photon- 
echo (PE) direction 57|. The frequency components of the pulses lie within the optical regime, so 
they can induce the transitions enumerated in the previous section. Traditionally, in the MDOS 
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literature, the intervals between the time centers of the pulses are called coherence t = t2 — ti, 
waiting T = t^ — t2, and echo t = t4 — t^ times, respectively. Here, ^4 is the time of detection of 
the signal 58|. We shall only consider rephasing photon-echo signals, where ti < t2 < ^3 < ^4, 



where the inhomogeneous broadening is rephased [59|. Due to these explicit interval dependences, 
the collected signal can be expressed as PpE{T,T,t). As explained in our previous study [34^, we 
may regard the PE experiment as a QPT of the single-exciton manifold dynamics of the dimer 
as a function of T. In fact, the polarization signal may be expressed as a linear combination of 
elements of the process matrix x{T)'- 



[PpE]el',e2',eleii.'t) " XI '^'^2 '^'^3-^61^62,63,64 W ' (13) 



p,q,r 



where, 



61,62,63,64 v"/ 

-{^^pg■ (ii)if^qa- e2)Ggp{r) 

^{[i^^ag ■ e3){^^ag ■ ei)Gag{t) 
^ iXggqp(T) — 6pq — XaaqpiT)) 

+ if^fi3 ■ e3)(M//3 ■ e4)gf^{t)xmpiT) 
+ ((M//3-e3)(Ai/«-e4)^j„(t) 

-(Mag ■ e3)(M/3g " ^4)^ ^g{t))Xpaqp{T)]} , (14) 



and, 



61,62,63,64 V''/ 
= -if^pg ■ ei){^igg ■ e2)ggp{t) 

x{[(M;39-e3)(M/39-e4)^/3g(t) 
^iXggqp{T) - 6pq - Xmp{T)) 
Ht^fa ■ e3)(M/a ■ (ii)Qfait)Xaaqp(T) 
Hi^^fa ■ e3)(M/,3 ■e4)^/;3(t) 

-(M/3g • e^){^^ag " ei)gc,g{t))Xal3qp{T)]} . (15) 
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The coefficients C^. for p G {a, /?} are frequency amplitudes of the laser pulse which is centered 
at Ui, evaluated at the transition energy Upg-. 

CI = --v/2^e-"'("-''-"^)'/^ (16) 
i 

and 



^..(r) = e(r)e(-*'^'^-^»^)" (17) 

is the propagator of the optical coherences in the coherence and echo times, which, has been 

taken to be the product of a coherent oscillatory term beating at a frequency Uij and an exponential 
decay with dephasing rate Tij. This propagator is defined only for r > via the step function 0(t). 
The frequencies of the coherences in the coherence and echo intervals have opposite signs, reflecting 
the rephasing character of the signal. In optical PE experiments, it is customary to assume that 
the free- induct ion decay characterized by the evolution of optical coherences in the coherence and 
echo times is well characterized, and given by expressions of the form ( ITTl) . The reason is that 
the characteristic energetic scales of the vibrational degrees of freedom are much lower than the 
optical gap, so the only nonunitary dynamics they induce in the optical coherence is, to a good 
appr oximation, restricted to pure dephasing which can be inferred from the polarization signal 



103] . The dynamics in the waiting time is more complex, consisting of small frequencies due to 
excitonic superpositions which are strongly influenced by the bath. It is the latter interval where 
QPT will prove most useful. 

The polarization signal yields a linear combination of elements Xabcd{T) weighted by the prob- 
ability amplitude to prepare a state \c){d\ with the ffist two pulses and detect \a){b\ with the 
third pulse and the fourth heterodyning pulse. These probability amplitudes can be controlled 
by manipulating the polarization of the pulses and the frequency amplitudes for the resonant 
transitions C^^, C'^ , C^^. In essence, state preparation and QST are implicit in the coherence and 



echo times (see 



34|). In a different context, Gelin and Kosov had previously hinted at a similar 



idea by identifying these times as "doorway" and "window" intervals 60|. By conducting several 
experiments varying these control knobs and collecting the signal from each of these settings, a 
system of linear equations can be established whereby the elements of x{T) can be inverted, and 
therefore QPT is achieved. This statement is correct provided that besides the free-induction decay 
rates Fjj, the parameters ujag-, UJjSg, /-i^g, /-i^g, fi^g, and fi^g are all known or can be obtained self- 
consistently during the experiment. We will elaborate on these points for the case of a homodimer 
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in section V. 

Notice that Eqs. (fT^ . (fT^ . and (fT^ monitor all the 12 real valued paramenters involving 
Xabcd{T) for a, b,c,d ^ {a, P}, so that they allow for the QPT of the single exciton manifold, which 
is an effective quantum bit (qubit) system. However, we have also kept track of the elements 
Xggcd{T) c,d & {«,/?}, that is, the possibility of amplitude leakage errors from the single-exciton 
dynamics channel to \g) {g\. It is known that whereas the excitonic dynamics occurs in femtosecond 
timescales, exciton recombination happens in the order of nanoseconds. Therefore, these decay 
channels could be potentially ignored in many experimental systems. We shall keep them in our 
theoretical analysis as they do not increase the complexity of the problem by much, although in 
situations where this could be problematic, we could accordingly disregard them. 



IV. QPT FROM 2D SPECTRUM OF PE 

As mentioned, QPT can be carried out from data resulting from a series of experiments varying 
colors and polarizations of the pulses. The necessary information can in principle be obtained 
by collecting a single point for a fixed pair of r and t points for each of the experiments. Often, 
however, the PE signal is collected across many r, T, t points, and conveniently processed into a 
2D correlation spectrum in the conjugate frequency variables Ur and uf. 



S{Ur,T,Ut) 



which still evolves in the T coordinate 



nj — I Jo 



;i8) 



104| . By performing the integrals of Eq. (fTSjl using Eqs 



USD, ([H]), and ([ISD, we obtain: 



S{uJr,T,UJt) = i l^^rnMh,ni^t)SmniT). (19) 

m,n=o,/3 

The spectrum consists of a sum of four resonances at {uJr,uJt) G 
{{uag,<^ag), {^ag,'^l3g), {^I3g,<^ag), {^I3g,'^l3g)}, which correspoud to the frequencies of the op- 
tical coherences at the coherence and echo times. These resonances are modulated by lineshape 
functions of the form. 



lr,mM = ■( _ _-r V (20) 

lyOJ-p ^mg ^J- mg ) 

kni^t) = (21) 

l[-UJt + Ung - U ng) 



11 



105 |. 



that correspond to the one-sided Fourier transform of the propagator along each r and t axis 
The peaks are centered about u = Umg and have a width parameter Tmg- The difference in signs for 
the Fourier transform in Eq. f |T8l) guarantees that all the resonances appear in the first quadrant of 
both frequency axes. The expressions for the amplitudes Smn{T), associated with peaks centered 



at {uj-r,ujt) = {uJmg,(^ng), are givcu by 



106| : 



^{Cuslit^ag ■ 63) (Mag " ^AjiXggaaiT) - 1 

+C<^3[(^/" ■ e3)(M//3 • 64) - (M/3<; ■ e3)(Mag " 64) (T)] } 

Xaaf}a 

+ ■ e3)(M//3 ■ e4)X,3/3/3a(T)] 

+Cl^M^^f<^ ■ e3)(M//3 • 64) - {fJ-pg ■ e^){^i^g ■ ei))xappa{T)]} (22) 



S^p{T) = -^C,Q(M.g-ei)(M.g-e2) 

^{Cmf^Pg ■ e3){flf^g ■ e^){Xggaa{T) - 1 - XfS^aaiT)) 

+ (M/a •e3)(M/ in 

+CZ.M^^fP ■ e3)(M/a ■ 64) - {flag " e3)(M/3g " 64) )X^aaa (T)] } 
x{C2J(A^/35 ■ e3)(M/3g ■ e4){X99^a{T) - X/3/3/3a(T)) 

+ (M/a-e3)(M/ in 

+C:^M^^fP ■ e3)(M/a ■ 64) - (Mag ■ e3)(M/3g ■ e4))X/3a/3a(T)]} (23) 
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x{Cf3[(/^/39 ■ es) (At/33 ■ e4){x99Pp{T) - 1 - xpm{T)) 
-^<Q(M/33-ei)(/x,3-e2) 

in 

+Cm{t^ffs ■ e3)(M/„ ■ 64) - (/X,, ■ es){fipg ■ e4))x/3aa/3(T)]} (24) 



= -^Cf,C7f^(^;33-ei)(M/33-e2) 

x{C"3[(Mag ■ e3)(A^ag ■ e4)(X99/3/3(^) - 1 - Xaal3piT)) 

+ ■ e3)(Ai//3 • ei)x(3i3f}^iT)] 

+Cl^M^^fc^ ■ e3)(M//3 ■ 64) - (M/3g ■ e3){flag " 64) )Xa/3/3/3 (T)] } 

-^C7^iQW.-ei)(M.3-e2) 

x{C"3[(Ma(, ■ 63) (Mag ■ e4)(Xgga/3(T) (T)) 

Hl^fP ■ e3)(M//3 ■ ei)xppap{T)] 

+C^M^^fo. ■ e3)(M//3 • 64) - (A^/33 ■ e3)(Ma3 ' f^i))XaPap{T)]} (25) 

Typically, the probed samples are in solution, so the molecules in the ensemble are isotropically 
distributed. The isotropic average (■) for a tetradic (/x„ ■ ei)(/x^ ■ 62) (^t^ ■ 63) (/x^ ■ 64) is given by 



((/^a ■ ei)(/x;, • e2)(Aie " e3)(Aid • 64)) 

/ J eie2e3e4;mim2m3m4 
mim2»n3m4 

X[(A^a ■ ^l)(Mfe • »T^2)(Mc ■ "^3)(Md ■ "^4)], 



(26) 
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30 



ei 626364 ;mim2m3m4 
1 

TT7[ ['^6162 "^6364 1^6163 '^6264 



'^ei64'^( 



62 63 J 



X 



^mim2^m,3'm4 



^mi m3 014 
^mi 7714 ^m2'm3 



(27) 



4 -1 -1 
-1 4 -1 
-1 -1 4 

where Cj and mj are the polarizations of the pulses in the lab and the molecular frame, respectively. 
The isotropic average consists of a sum of molecular frame products [(/x^-mi)(^i{,-m2)(/^c-Tn3)(/Xj^- 
m.4)] weighted by the isotropically invariant tensor /ife26364;mim2m3m4. 

Since the information in Eqs. (1131) . (1141) . (fT5|) is in principle contained in Eqs. (fT5|) . (IT^ . and 
([T^ . several conclusions from our previous study are immediately transferable: The elements of 
x{T) contained in Eqs. (122|) . (|23|) . (IMj) . and ( 125|) can be all be extracted by repeating a number of 
experiments with different polarization configurations for the fields and two different waveforms for 
the pulses. Under different motivations, theoretical proposals for manipulating 2D-ES using pulse- 
shaping capabilities have been previously reported 62|, |63|]. An extensive study of this possibility 
for a heterodimer will be presented in the second article accompanying this study. 

Eqs. (I22p . (125]) . (I24p . and ( 12 5 p can also be derived by book-keeping the double-sided Feynman 
diagrams that oscillate at the particular frequencies for the coherence and waiting times in each of 
the four resonances (we refer the reader to Fig. 1). In analyzing the possible pathways in Liouville 
space, we make use of the rotating-wave approximation (RWA): Perturbations which are propor- 
tional to e^**'»'^+*'^'(*~**)/i-ej can excite the ket and de-excite the bra, whereas the ones proportional 
to e^^i''^-'^^i^^^^i) jx - Si can deexcite the ket and excite the bra. As an illustration, consider the signal 
SajsiT), which arises from diagrams oscillating with frequency Uga at the coherence time and 
at the echo time (Fig. 1(b)). The two states at the coherence time which can oscillate at Uga 
are or but the latter cannot be produced by a single action of the dipole operator 

on the initial ground state \g){g\. Therefore, |(7)(a;| is the only possible state for the coherence 
interval, and is produced by acting the first pulse on the bra of the ground state: \g){g\ — ?■ 
Similar considerations imply that the state at the echo time must be \(3){g\ or Given these 

constraints, we are ready to enumerate the possible initial and final states for the waiting time 
interval which are compatible with these restrictions. By exciting the ket or deexciting the bra 
of with the second pulse, the following initial states \c){d\ for the quantum channel can be 

produced: \c){d\ G {|Q;)(a;|, |/3)(q;|, The final states \a){b\ G {|Q;)(a|, |/3)(q;|, |5')(5f|} 
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Figure 1: Liouville space pathways corresponding to each of the four resonances in the 2D-ES of a coupled 
dimer. The amphtude Smn{T) corresponds to the peak located at {(jOt^ijOi) = (urng^^ng), which are the 
values of the optical frequencies at the coherence and echo time intervals r and t, respectively. These 
amplitudes provide information on the coherent and incoherent excitonic processes at the waiting time 
T by enumerating the possible initial and final states at the waiting time T which satisfy the PE phase 
matching condition for the pulses acting at times ti, t2) ^3) ^4- The information contained in the amplitudes 
Smn(T) can be distilled to reconstruct the process matrix for the single-exciton manifold of the dimer, 
thus allowing for a QPT. 

can all give rise to \P){g\ or |/)(a;| by exciting the ket or deexciting the bra with the third pulse. 
Therefore, in principle, there are 4 x 3 = 12 possibilities for Xabcd{T) which can be detected in at 
SajsiT). However, we assume that the state \g){g\ does not evolve to other states due to the bath: 



XabggiT) = SagStg, (28) 

This assumption is quite reasonable, as we are ignoring processes where phonons can induce optical 
excitations from \g){g\. This condition is present in Eqs. 0141) . and f[T51) in the 5— function terms 
and in Eqs. (122|) . fl23|) . flM|) . and f l25|) in the "-1" terms, which correspond to —Xgggg{T). This 
leaves 12 — 4 = 8 possibilities for Xabcd{T)- 



To be more explicit, consider the pathways in Sap{T) that monitor the population to coherence 
process Xi3aaa{T). Thcsc are displayed in Fig. (2). The pathway on the left involves represents an 
excited state absorption (ESA) from the single-exciton manifold, and is proportional to {—C"_^fi^g- 
ei)(C^2Mag ■ ^2){C!^^fJ-fi3 " ^z){t^fa ' ^4); cxprcssiou which can be immediately read out from the 
diagram: Each interaction with the field picks up a factor corresponding to the amplitude of the 
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Figure 2: A more detailed view on the Liouville space pathways corresponding to the monitoring of 
the population to coherence process Xl3aaa{T) in the peak at [uj^-^i^t) = {'^ag^'^pg) of the 2D-ES. These 
diagrams belong to the amplitude Sap{T) and can be easy constructed by taking into account the PE 
phase-matching and the resonant conditions. 

transition, which depends on the ahgnment of the corresponding dipole with the polarization of 
the pulse, as well as the frequency amplitude of the pulse at the given transition. A minus sign is 
included if the perturbation is on the bra. Similarly, the pathway on the right involves stimulated 

emission (SE) and is proportional to (-C^^Mag " ei){CZ^ii^g ■ e2){-C^^lJiag ■ f^^){t^pg ■ 64)- 

The rest of the diagrams for all the peaks can be systematically analyzed in the way described 
above. In general, the pathways we need to consider can be classified in ESA, SE, and ground state 
bleaching (GSB) processes. GSB processes are the ones that take \g){g\ at the end of the waiting 
time to a dipole active coherence involving an excited state. ESA pathways, which are proportional 
to dipole transitions involving the excited state, differ in sign from SE and GSB pathways, as can 
be easily seen by inspection. 



Fig. (3) provides a mnemonic device to keep track of the Liouville pathways that each peak in 
the 2D electronic spectrum monitors, and therefore, also provides a scheme of the QPT protocol. 
The ujr axis can be associated with a particular state preparation whereas the Ut axis with a 
particular detection. Each peak reflects a nontrivial number of processes in Liouville space. As 
an illustration (see Fig. (4)), we consider the ideal case where the bath does not interact with 
the system, in which case, a very simple picture is recovered: The off-diagonal peaks beat at the 
coherence frequency and the diagonals remain static. This case can be easily derived from Eqs. 
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Figure 3: Summary of QPT for a coupled dimer in the 2D-ES. The Liouville pathways depicted in Fig. 1 
can be condensed into this diagram. The horizontal axis for the coherence frequency uJr is associated with a 
state preparation, whereas the vertical axis for the echo frequency ujt corresponds to a detection. The four 



resonances labeled as {m,n) correspond to peaks located at {00-^,001) 



[ijj. 



mgi 



^ng)- Their amplitudes con- 



tain information on Xabcd(T), where cd is the state prepared at the beginning of the waiting time interval, 
and ab the state detected at the end of the sa,in6 interval. For instance, the peak at (wct, uiis) keeps track 
of the elements Xabcd{T) where \c){d\ G {|q)(q|, |/3)(q|, |5)(c/|} and \a){b\ G {\/3){a\,\f3){f3\,\a){a\,\g){g\}. 

fl22|) . fl23|) . fl2i|) . and (ES]) by substituting Xabcd{T) = dabbed + ^ac^Me"''^""^, that is, populations 
remain static whereas coherences beat at difference frequencies. 



V. THE CASE OF THE HOMODIMER 

To gain insights into the described QPT protocol, we specialize the results above to a coupled 
homodimer. In the following subsections, we discuss, for this particular case, (A) the Hamiltonian 
and the transition dipole moments involved in the experiments, (B) properties of the spectroscopic 
signals under isotropic averaging, (C) stability of the numerical inversion, (D) analytical expressions 
of the elements of x{T) in terms of the peak amplitudes of the spectra, (E) a procedure to extract 
the angle between the dipoles, (F) a summary of the QPT procedure, and (G) a numerical 
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Figure 4: 2D-ES of a coupled dimer in the absence of interactions with a bath. Under unitary dynamics 
of the excitonic system, each of the four resonances keep track of the elements of x{T) indicated in the 
diagram. Notice that diagonal peaks do not oscillate as a function of waiting time T, whereas off-diagonals 
beat at the frequency equal to the difference in energies of the single exciton eigenstates. 

example with a model system. A similar study focused on the heterodimer will be presented in 
the second part of this series. 

A. Hamiltonian and transition dipole moments 

In the homodimer, the two sites are identical chromophores with energies ua = ^^b = ^, and 
the Hamiltonian in Eq. ([5]) and ([7]) is given by: 



which we have diagonalized with the symmetric a'^\g) and antisymmetric a'^\g) single exciton 
states given by: 




(29) 
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1 

71' 

1 



4 = (30) 

The splitting between the two delocahzed states is just 2 J. Using Eq. (fTOj) . the transition dipoles 
take the simple forms: 



t^ag — 








^'P9 = 




^/a = 


^ag^ 


= 





A + ds) 
A - ds) 



(31) 

Interestingly, these expressions are independent of the coupling J. Also, notice that /x^,^ and 
are perpendicular to ^i^g and ^if^ (see Fig. 5). Denoting the norm of each site dipole by 



\dA\ = \dB\=d, (32) 



the following relationships follow: 



\^^ag\ = IM/oI = fl'ag = V2(i COS 

= IM//?! = f^i3g = V2dsm . (33) 

As expected, in the degenerate limit that = or vr (the site dipoles are parallel or antiparallel) , 
one of the delocalized excitons becomes dark and there is only one bright transition from the 
ground state. If this is not the case, in general, the two transitions are bright and their dipoles 
perpendicular to each other. Furthermore, as a difference with the heterodimer case, there are 
only three (instead of four) different transition dipoles in the homodimer, and two of them are just 
negatives of each other. The degenerate case will be discussed as a limit of the more general one 
in the next paragraphs. 
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(a) (b) 




Figure 5: Transition dipole moments of a homodimer. Diagrams for (a) sites and (b) eigenstates. 
B. Isotropic averaging of signals 

An important observation regarding isotropic averaging follows: 

Claim.- Upon isotropic averaging, signals stemming from coherence to population or popula- 
tion to coherence transfer cannot be monitored in the 2D-PE spectrum of a homodimer. 

Proof.- For simplicity, align /x^^ and in the y and z directions in the frame of the molecule 
and use Eq. ( l26l) to argue that the only terms in the sum that contribute to an isotropic averaging 
are the ones where only two or four polarizations of the field are the same. In Eqs. fl22l) . fl23|) . 

and (125|) . all the dipole-polarization terms corresponding to coherence to population and the 
opposite processes involve three dipoles of the same kind and a perpendicular one. Therefore, they 
vanish under isotropic averaging. As an example, consider the the terms associated with Xi3aaa{T) 
in Sai3(T): 
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= 0, 

((-Q/^ag ■ ei)(C2/^ag " e2)(-C"3Mag " e3)(M/3g " 64)) 
« ((Mag • ei)(Mag " 62) (^ag " e3)(M/3g " 64)) 

= 0. 

□ 

The claim above allows for a considerable simplification of Eqs. (12^ . fl25]) . ([21]), and 025 p . Each 
of the peaks in the 2D spectrum keeps track of less elements of the process matrix x{T) upon 
isotropic averaging: Only population to population and coherence to coherence transfers can be 
monitored. The results for the peaks of an isotropically averaged 2D spectra are presented below. 
We have taken the shortcut notation (■)eie2e3e4 5 which denotes the isotropically averaged signal 
stemming from the two pulse polarization configurations (ei, 62, 63, 64) = {z, z, z, z), (2, z, x, x), 
so that the terms (5'mn(r))eie2e3e4 and {S{uT-,T,Ut))eie2e3e4, have the obvious meanings. 
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TABLE 1. Isotropically averaged 2D-ES peak amplitudes for the zzzz configuration 



('S'a/3 (2^)) zzzz 



ZZZZ 

^llf^kiXgg^eiT) - 1 - X/3/3/3/3(T)) 
-iCt,CZ^CZ,{{-)^txiy^,X?o.a0{T)] 



X [^A*ig(Xffffaa(r) - 1 - Xaaaa(I')) 
+ ^MasM|gX/3/3ac«(r)] 

-iczAcL\i.-)^i^lAgx<^f^?<^{n 



= -iC^ C^ C" 



^[Y^f4gt^lg(Xggfll3{T) - 1 - Xaocl3/3{T)) 



-iCt.CSAii-)^4gt'lgX<^^MT)] 
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TABLE 2. Isotropically averaged 2D-ES peak amplitudes for the zzxx configuration 



zzxx 
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><[-^f^lgf^%{Xggaa{T) - 1 - X^^aa{T)) 
(T)] 



x{Y^^^h(xggf^p{T) - 1 - xpmiT)) 

2 

CS, C^s [(+) 3^M^,/^a,X/3aa/3 (T)] 



•/oa /^a 



^[-11- f^ag iXggaa {T) 1 Xceaaa (^) ) 



'15' 
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{Spain, 



A^^^Pg^^agiXggm{T) - 1 - XaaPp{T)) 



+^^4gXf}miT)\ 
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We focus our attention on experiments with short pulses that are broadband enough to create 
either exciton |a) or with the same amphtude, that is, C^^ = C, for a purely imaginary constant 
C, for all p and Ui. This condition can be easily relaxed, but we shall proceed with it to analyze 
our QPT protocol with more detail. Using the condition of Eq. ([3]), we can eliminate the variables 
XggaaiT) and XggfspiT) for Xaaaa{T) , Xi3f}aaiT) , Xaai3i3iT) , Xi3f3i3i3{T) ■ Also, taking advantage of Eq. 
(II]), we discard XpapaiT) and Xcy-PPaiT) and keep Xa/SaisiT) and XpaapiT). From the left hand side 
(LHS) of the real and imaginary parts of the spectrum (see Tables 1 and 2), we derive the following 
real valued matrix equation: 



1 .A 



1 ,,2 



5 Mc 
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Similarly, the right hand side (RHS) of the spectra yields: 



Xcccc(T) 
«{X3oo,3(T)} 



3*{<Scc(T))„„} 
K{(S„„(T)>„,,} 

9{<Scc.(T))„„} 
a{<S„„(T)>„,,} 

a{(s„;3m>„„} 



(31) 
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Inverting Eqs. and ( l35l) yields most of the elements of x(^) involving the single-exciton 

manifold. Whereas the presented QPT for the homodimer is partial, no complicated pulse shaping 
efforts need to be carried out. Instead, the requirement is standard pulse polarization control 
achievable with current experimental capabilities fill. [12I. 64-69|. 



The transition dipole moments must be well characterized in order to construct the matrices in 
Eqs. (l34|l and ( !35|) . This requirement is self-consistently fulfilled by collecting the spectra in the 
collinear and cross-polarized configurations. Notice that Xaaaa{T) and Xpi3aa{T) are exclusively 
monitored in the LHS, and Xi3i3i3i3{T) and Xaai3i3{T) only detected in the RHS. However, coherence 
transfer terms Xai3ai3{T) and Xi3oiap{T) are repeatedly monitored in different peaks in both sides of 
the spectra. Due to this repetition, there are redudant equations that allow for the self-consistent 
extraction of the angle (j) without compromising the inversion of the elements of x{T)- Details 
about this parameter extraction are developed in subsection E. For the time being, we assume 
that the information about the transition dipoles is previously known. 
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C. Stability of the Quantum Process Tomography protocol for a homodimer 



In order to characterize the stability of inversion of x{T), we can arrange Eqs. f lMl) and f l35l) 
into a single matrix equation Mx{T) = S{T), where M is a 16 x 8 matrix of dipole moments, 
x(T) is a vector of 8 unknowns, and S{T) is a vector of 16 real valued am plit udes extracted from 
the signal. Denoting with || ■ || the spectral norm of a vector or a matrix [TOj, we obtain a bound 
on the relative error of the inverted vector x{T) which yields the QPT: 

l|Ax(r)|| ^ II AS{T) II 

where Ax{T) and AS{T) denote errors in x{T) and S(T) upon inversion, respectively, and the 
condition number k, is given by: 



K =11 M nil M"^ II (37) 

The lowest possible value for a condition number is k = 1. In Fig. (6), a plot of k vs. (p (red) 
is displayed. As expected, very large values of k, which denote unstable inversions, are expected 
for systems where the site dipoles are aligned or antialigned, when one of the eigenstates of Hs 
becomes dark. We also indicate a range of angles where k is below a reasonable threshold, say 
K < 15 (blue horizontal line), which consists of angles in the range O.Stt < < O.Tvr. Also, note 
that K is symmetric about the minimum of ft at = 7r/2, where the best inversion is carried out 
with K, ~ 3.9. 



D. Analytical expressions for x(^) 

To gain insights into the QPT protocol, we shall derive explicit expressions for the elements 
of x(T) in terms of the amplitudes of the spectra and the angle between the dipoles. First, 
we substitute Eqs. ( l32l) and (l33l) into the LHS of the spectra through Eq. ( IMl) . The following 
expressions are obtained after inverting the resulting matrix equation: 



(T) - 1 = '"'"^^/^^ 



(iC3)(20rf4) 

^'^{{Saa(T)) zzzz ~l~ '^{Saa(T^) zzxx 

+ {Sap{T))zzzz + '^{Sap{T))zzxx} (38) 
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Figure 6: Stability of the homodimer QPT protocol. As described in the text, the QPT protocol depends 
on the inversion of a matrix which is a function of transition dipole moments. The condition number of 
the matrix (k) vs. the angle between site dipoles (0) is plotted in red. The blue line plots the constant 
value of K = 15, as a reference to indicate that for the range of angles O.Svr < </> < O.Tvr, the condition 
number is below that value. 



(T) 



Xpi3aa{T) 



sec^(0/2) 
~(zC3)(20rf4) 

X3?{(COS0 - l){{Saa{T))zzzz + '^^{S aa{T)) zzxx 

+ COS (t){{S^p{T))zzzz + '2{S^p{T))zz..)} (39) 



sec^(0/2) 
(iC3)(20#) 

X^{C0S (l){{Saa{T))zzzz + '^{Sac,{T))zzxx) 

+ (cos0 + l){{S^p{T))zzzz + 2(S,;3(T)),,,,)} (40) 
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(iC3)(20#) 

+ (cot2 ^ + 2tan2 ^ + 5)(^«/3(T)),,,, 
-(3cot^ ^ + tan^ ^){Sai3(T)) zzxx} 



1 



(iC3)(20rf4 



x^{{cot^ I + 2tan2 | + 1)(5„„(T)), 
_(3 cot^ ^ + tan^ ^ + 8)(5„„(T)),,,, 

—2{{S al3{T)) zzzz + 2(S'Q^(T))22a;a;)} 



-^^{(csc^ ^ sec^ ^)(5'„^(T))2^a;a;} 



(zC3)(4rf4) 2 2 



^>{X/3aa/3(T)} = (^(g3^)(g^4) ^{(cSC^ ^ SeC^ ^) (^aa (T) ) 

^^{(csc^ ^ sec^ ^)(S'a„(r))^^^^} 



(iC3)(4#) 2 2 

Similarly, the RHS of the spectra through Eq. (l35l) yield: 



(zC3)(20#) 
x3?{(5^;3(T)),,,, + 2(^^^(T)),,, 



_ csc2(0/2) 



(zC3)(20#) 
x3fj{(cos0 + 1)((5^^(T)),,,, + 2(5^/3(T)), 



cos 



0((^/3a(r)),,,, + 2(S'^„(T)),,,,)} 
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(iC3)(20rf4) 
x3ft{cos0((5;3/3(T)),,,, + 2(5^;3(T)),,,,) 

+ (cos0 - + 2(V(^))— )} (47) 



x$R{2((^;3;3(T)),,,, + 2(5^^(T)),,,.,.) 
+ (tan2 ^ + 2cot2 ^ + 5)(5/3a(T)),,,, 

-(3tan2 ^ + cot^ ^)(5^„(T)),,,J (48) 



x3?{(tan2 ^ + 2cot2 ^ + 1)(5;3/3(T)),,,, 
-(Staii^ ^ + cot^ ^ + 8)(5^^(T)),,,, 

-2((V(T)),,,, + 2 (5;3a(T) ),,,,)} (49) 



^{Xa^apiT)} = - ^^^3^^^g^4^ ^{(seC^ |cSC^ |) (V (^) )^^^ J 

-5>{(sec^^csc^^)(S'/3<:,(T))2^^^} (50) 



(^C3)(4rf4) 2 2- 



Q'{X/3aa/3(T)} 



As expected, the RHS equations above can be obtained from the LHS equations upon the substi- 
tutions a — )■ /3, /3 — )■ a, and — i- vr — (see Eqs. (133|) ). 

From Eqs. ( l43|) and ( l50l) . we notice that the imaginary parts of 

(S'a/3(T))^^^^, (S'a/3(T))22a;a;, {S^a{T)) {S^a{T))zzxx are all exclusively proportional to 



-^{(sec^ ^ csc^ ^)(5/3/3(T))^^^J 



(2C3)(8rf4) 2 2' 

1^ -53{(sec2 I csc^ '^){Spp(T))zzxx} (51) 



(iC3)(4c/4) 2 2- 
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^{Xai3ai3{T)} ■ The least-squares solution for '^{XafSafsiT)} using Eqs. and (135|) is sim- 

ply the average of these four values; the analogous conclusion holds for Q{xi3aai3{T)}. The real 
parts, ^{xai3ai3{T)} , ^{xi3aai3{T)} Satisfy less trivial relationships. Each of them appears twice in 
the equations, once in the RHS and another time in the LHS of the spectra. It is this redundancy 
in the spectral information what allows for the extraction of the angle between the dipoles. 



E. Determination of the angle between the two dipoles 

Equating the values of ^{Xai3ai3(T)} from the expressions in Eqs. (l4T]) and (HHj) yield the 
following quadratic equation in ^ = tan^ (|): 



^{[{Sl3aiT))zzzz — '^{Sl3a{T))zzxx — '^{S ap{T)) zzzz + {Sali{T)) zzxx]^^ 
+ ^{Spa{T))zzzz + 2{Si3i3(T)) zzzz + 4:{S iS/siT)) zzxx 
— 5{Sctj3(T)) zzzz ~ '2{Saa(T)) zzzz ~ 4:{Saa(T)) zzxx]^ 

+ ['2{Si3a(T)) zzzz — {Si3a(T))zzxx — {S afiiT)) zzzz + "^{3 afiiT)) zzxx\] 

= 0. (52) 
A similar expression can be found for ^{xi3aai3{T)} from Eqs. (H2|) and f H9|) : 



^{[2{Saa(T)) zzzz — {SaaiT)) zzxx — {S jipiT)) zzzz + 'i{S jifiiT)) zzxx]^ 
~\~\{S aaiT)) zzzz ~ ^{Saa(T)) zzxx ^ ^{SajsiT)) zzzz ~ ^{Sa/siT)) zzxx 
~{Sfip{T))zzzz + 8{Si3j3(T))zzxx + 2{S pa{T)) zzzz + '^{S pa{T)) zzxx\i 
+ [{S aa{T)) zzzz — '^{Saa(T)) zzxx — 2,{S pii{T)) zzzz + {S pii{T)) zzxx\] 

= 0. (53) 

The identities in Eqs. (1521) and (I53p are remarkable in the sense that they are satisfied at every 
waiting time T: they do not depend on short time coherent dynamics. By monitoring the peak 
amplitudes from the spectra arising from the two different polarization configurations, the angle (j) 
between the two site dipoles can be readily extracted using either expression. This determination 
is robust because it can be repeated for every value of T for which the signal has been collected. 
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F. Summary of Quantum Process Tomography protocol for a coupled homodimer 

We proceed to summarize the algorithm of the QPT protocol for a coupled homodimer. 

1. Obtain the amplitudes {Smn{T)) eie2ese4 for m,n & {a,/?} and the two polarization settings 
(61,62,63,64) = {z, z, z, z), {x,x, z, z). This information can be extracted from the two 
respective polarization controlled 2D-ES, {S{ur,T,ut))eie2e3e4- For simplicity, all the pulses 
are taken to be of the same duration (short compared to the timescales of excited state 
dynamics) . 

2. Extract the angle from the data from step 1 and Eqs. (l52l) and (l53l) . 

3. Plug in information obtained from step 1 and the angle (p from step 2 into the expressions for 
the elements of x{T) in Eqs. ( l38l) -fl5Tl). Some important observations: (a) These expressions 
are all proportional to the factor {C^d"^)'^. The norm of the dipole d can be extracted from 
the intensity of the absorption spectrum of the monomer. If this information is not readily 
available, the results are known up to this constant factor, (b) By construction from Eq. 
(152|) . the calculated value of ^{Xai3ai3{T)} will be the same using either Eqs. (HT]) or (H8|) . 
The same holds for Eq. ([53]), ^{xi3aap{T)} and Eqs. fHTj) and ipH]). 



G. Numerical example 

In this subsection, we illustrate the described QPT protocol with a model homodimer. Mar- 
cus and coworkers have recently reported a synthetic system of porphyrin molecules which self- 
assemble into homodimers under the presence of liposomes The parameters of this system, ex- 
tracted from phase-modulation electronic coherence spectroscopy (PM-ECS), are u = 16633 cm~^, 
J = 175 cm~^, and = 65°. The transition energies to the eigenstates are Uag = 16458 cm~^ 
and Ui3g = 16808 cm~^. Information on the spectral density of this system is not available in the 
literature yet. We adopt a simple system-bath model based on the secular Redfield approach and 
independent bath approximation (IBA) for each site. The weak system-bath model is reasonable 
since porphyrins are rigid molecules which change their structures minimally upon electronic ex- 
citation. The IBA must be reexamined, since the liposome media guarantee a bosonic bath that 
could be strongly correlated in both sites. Nevertheless, the purpose of this example is not to 
provide an exact account of the excited state dynamics of this system, but rather an illustration 
of the QPT protocol using reasonable timescales that one might encounter in a realistic setting. 
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A careful study of the precise bath-induced dynamics in this system is beyond the scope of this 
study and shall be addressed in future work in collaboration with an experimental realization. 

We consider a harmonic bath with an Ohmic spectral density: J (to) = -^uje^^l'^" ^ with 
A = 100 cm~^, = 150 cm~^ at a temperature T = 273 i^. Identical baths are assumed to 
be diagonally and linearly coupled to each of the sites. We closely follow the calculation reported 
in 72| and adapted in [34] . The dynamics of the total excitonic system, which is a proper density 



matrix, is governed by the following equation of motion: 



p{T) = -Aiis.piT)\-Tlp(T) (54) 

where 7^ denotes the time-independent sparse dissipative superoperator containing only a few 
non-zero elements listed in Table 3. Since p(T) only depends on p(T) and not on the value of the 
quantum state at previous times, the simulated dynamics are Markovian. 



TABLE 3. Values (in j s ^) of non-zero rates of the secular Redfield tensor 





8.02 X 10" 


-4 


R an = e~'^°'^^^B'^ Baa 


5.07 X 10- 


-3 


Ra/SalB = RpaPa 


2.93 X 10- 


-3 


Ragag Rgaga Rfafa Rafaf 


1.23 X 10" 


~2 


RpgPg = RgPgP = Rfl3fl3 = Rfiffif 


1.45 X 10" 


-2 


^fgfg = ^gfgf 


4.77 X 10" 


-2 



It is well known that the secular Redfield equations guarantee thermal equilibrium since the 
population transfer rates satisfy Raapp/ RpiSaa = e""'^"''/^-^"^, where ks is the Boltzmann constant. 
Also, Rfgfg will not be relevant for the calculations, as coherences between the ground state and the 
biexciton are never created in the PE experiment. The free-induction decay rates for the coherence 
and echo intervals will be taken for simplicity to be the same, T^g ~ f» \{Ragag + R/Sg/^g)- This 
restriction is by no means necessary, but will simplify the simulations below. 

The non-zero elements of x{T) for the single-exciton manifold are presented in Table 4. 
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TABLE 4. Nonzero elements of x{T) involving single-exciton states for the secular Redfield model 



XaaaaiT) 1 _ g-^/S/^-^ 



In this particular calculation, coupling to the photon bath has been ignored beyond the ul- 
trashort pulses, as spontaneous emission occurs in the order of nanoseconds, i.e. Xggaa{T) = 
Xggi3p{T) = 0. Therefore, Eq. ([3]) can be readily confirmed (T) = 

Xaa(3f3iT) + Xi3i3aa{T) = 1. Also, duc to the sccular approximation, X/SaafsiT) = (Xa/sf^aiT))* = 0. 
The picture of the secular Redfield equations is very simple and provides transparent means for 
understanding the QPT protocol for the homodimer: The evolution of populations and the co- 
herences independently satisfy standard first-order kinetic equations, leading to multiexponential 
integrated dynamics. 

In Fig. 7, we display the calculated 2D-ES of this model system. We consider the three pulses 
to be identical, centered about Ui = UJ2 = 003 = 16546 cm~^, of FWHM = 20 fs, i.e., a = 8.49 fs, 
which amount to an equal excitation amplitude C for both |a) and Rabi oscillations for a 
coherent superposition between \a) and occur with a period = 47.5 /s. We present several 
snapshots of the real and imaginary parts of the spectra at values of waiting time T corresponding 
to multiples of Tc/2, skipping Tc = 0, as our theory has avoided pulse overlap effects. In principle, 
as Eqs. (123|) and fl25|) indicate, the excitonic quantum beats associated with the term Xai3ai3{T) 
can be monitored by looking at either cross peak of the spectra. This feature is subtly manifested 
in every column of the figure, but more easily perceived in the real part the zzxx spectrum, where 
the peak at {uji3,uja) changes from red to yellow/green every interval Tc/2 before the bosonic bath 
has washed out significant portion of the coherent dynamics at about T = 5Tc. Also, incoherent 
population transfer primarily from to \a) (downhill) manifests as a decrease in amplitude of 
the peak at (upyUp) due to ESA and an increase in {uJi3,Ua) due to SE [2]. This effect can also be 
more obviously seen in the real part of zzxx spectrum. 



From the simulated spectra {S{uJr,T,Ut))eie2e3e4,, the extraction of the terms {Smn{T)) eie2e3e4, 
is achieved with high fidelity (>99%) by a nonlinear optimization routine based on the simplex 
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(a) Tc/2 



(b)Tc 



(c) 3Tc/2 



(d) 2Tc 



(e) 9Tc/2 



{f)5Tc 




1,62 1.64 1,66 1 



,62 1.64 1,66 1.68 1,7 



Figure 7: 2D-ES for coupled porphyrin homodimer with secular Redfield model. From left to right, 
we show the real and imaginary parts of the spectrum with zzzz polarization setting (first and second 
columns), and with zzxx setting (third and fourth columns). Each row represents a particular waiting 
time T, corresponding to (a) Tc/2, (b) T^, (c) 3Tc/2, (d) 2Tc, (e) 9Tc/2, (f) 5Tc, where = 47.5 /s is 
the period for one Rabi oscillation between \a) and |/3). The colormap is such that red is associated with 
positive numbers, green with values about zero, and blue with negative numbers. 
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Figure 8: Traces of {Smn{T)) e^eie^ei for m,n G {a,/3} and (61,62,63,64) G {{z,z,z,z), {z,z,x,x)}. (a) 
Real and (b) imaginary parts for the zzzz configuration; (c) real and (d) imaginary parts for the zzxx 
configuration. Each of the plots shows the evolution of the peak amplitudes SaaiT) (red crosses), Sa/siT) 
(green circles), SjSaiT) (blue squares), and Spp{T) (black diamonds). 



search method with bound constraints 73|. The signals are fitted to a sum of four different 
resonances as in the isotropically averaged version of Eq. ( IT9|) . The parameters Umg-, ^^ng, T, 
and (S'mn(T)) £1626364 are reconstructed from 2D-ES with a grid spacing of Aur = Aut = lcm~^ 
and a grid size of 1050 cm~^ for every axis. We present the results of this calculation in Figure 
8. Notice that the imaginary parts of the diagonal peaks are zero since no terms of the form 
XpaajsiT) = {xai3i3a(T))*a.Te Considered in the secular Redfield theory, and population transfer 
terms are purely real. 



Eq. fl52l) is solved at every waiting time T G [Tc/2, lOTc] yielding the roots ^ = 1.000, 0.4059 for 
every T, without variance after the fourth decimal digit, indicating its robustness for the inversion 
of (j). The same exercise with Eq. ( l53l) gives ^ = —0.4059, 0.4059. These values of ^ imply 
= 90°, 65°, 65°, 86.3i°, respectively. The last value can be discarded for it is not even a real 
number. The value of must be a root of both equations. Therefore, we can also discard 90°, 
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Figure 9: Extractable elements of x(^) for homodimer. (a) Processes starting from |a)(a|: Xggaa{T) (red 

crosses), \aa.aa 

(T) (green circles), Xi3i3aa{T) (blue squares), (b) Processes starting from \/3){/3\: XggHpiT) 
(red crosses), XaappiT) (green circles), Xpppp (blue squares), (c) Processes indicating coherence trans- 
fer: 3?{xa/3a/3(T)} (red crosses), '^{Xa.pa.piT)} (green circles), lk{xpaai3{T)} (blue squares), '^{xpaapiT)} 
(black diamonds). 

since it is only solution of the first equation, but not of the second. The result = 65° follows 
unambiguously, as expected. 

Finally, the terms (S'^.n (7")) £1626364 and the angle allow for the evaluation of the elements 
of x{T) which are extractable for the homodimer. Figure 9 shows that this reconstruction 
coincides with the analytical expressions presented in Table 4. The population decay terms 
Xaaaa{T),Xi3i3i3p{T) both start at 1 and reach exponentially, the second faster than the first, 
since is the excitonic state of higher energy. The population transfer terms XaappiT)-, XppaaiT) 
are complementary to the former ones, with the transfer from to \a) being faster for the same 
reasons just mentioned. The coherence term decays exponentially, with real and imaginary parts 
7r/2 phase shifted one from another. The calculated timescale of this decay (hundreds of femtosec- 
onds) is similar to the one inferred from the experiment reported by Lee and coworkers, where a 
superposition of excitons in the bacteriopheophytin and bacteriochlorophyll sites in the reaction 



center of purple bacteria is monitored indirectly through a two-color experiment 



74|. 
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VI. DISCUSSION 



In the present article, we have outhned a general theory for carrying out a QPT for a molecular 
dimer using the information contained in various frequency and polarization controlled 2D-ES. We 
started by providing the basic concepts of QPT, and operationally defined a QPT as a protocol 
to extract the process matrix x{T), which in principle completely characterizes a quantum black 
box, in our case, the box being the single-exciton manifold of the dimer. After reviewing the 
model Hamiltonian as well as the transition dipole moments of an excitonic dimer, we adapted 
the QPT theory presented in our previous work, where the nonlinear polarization was analyzed 
in real time for single values of r and t times (see Eqs. f[T^ . (IT^ . and (IT^ ) 34| . to the more 
standard and visual Fourier transformed signal collected along several values of these interval 
times. The central result of this exercise was Eqs. (122|) . (123|) . (^^, and ( !25|) . which from a purist 
stanpoint completes the QPT effort: The peaks in a heterodyne-detected 2D-ES can be expressed 
as linear combinations of elements of the process matrix x{T)- This information can be distilled by 
carrying out several experiments alternating the frequency components of the pulses as well as their 
polarization. By setting up a system of linear equations with this data, a linear algebraic routine 
yields the inversion of x{T) for every waiting time T. In order to get a more intuitive picture 
of this procedure, the particular case of a homodimer was studied in detail. The degeneracies of 
this system yield a perpendicular set of transition dipole moments which considerably simplify 
the theory (see Fig. 5). It was shown that under isotropic average of the signal, no population 
to coherence processes or viceversa can be monitored, impeding a full QPT for the single-exciton 
manifold of this system. However, the partially achievable QPT is very simple, robust with respect 
to transition dipole moment parameters as long as they are not aligned or antialigned, and readily 
implemented without pulse shaping. The only requirement is the collection of two polarization 
controlled 2D-ES. Numerical examples with a model homodimer validated the presented theory. 

The possibilities that QPT opens for the study of excited state dynamics in condensed molecular 
systems are as vast as the information acquired at the amplitude level of the evolving quantum 
state of the probed system. On the one hand, with the peaks in the 2D spectrum indicating a 
plethora of pathways in Liouville space, understanding of the dynamics is undoubtedly enhanced 
by the dissection of these peaks into processes described by the x{T) matrix. On the other hand, 
a plethora of questions can be addressed with this information, for instance: Is a Markovian 



36 



description accurate? 



79 



75 



78| If not, what is the degree of non-Markovianity of the dynamics? 



80| If it is Markovian, is the secular Redfield equation appropriate or are non-secular processes 

enient in the quantum states produced in the single 
83| What is the rate of decoherence of a quantum 



important? [81]. Is there any degree of entan 
exciton manifold upon photoexcitation? [82 



superposition between excitonic states? j84j-|88 



A few aspects have not been fully addressed with respect to the implementation of QPT of 
a molecular dimer. These issues will be carefully studied in future publications in collaboration 
with experimental groups. The role of static disorder in the eigenenergies of the system as well 
as in the distribution of the angle will necessarily yield an inhomogenously averaged signal 
from which the relevant information must be carefully extracted. We anticipate this feature to 
add another step of parameter fitting, but not change the results of our theory dramatically. 
Furthermore, we have ignored the possibility of resolving the vibronic structure accompanying 
each of the four resonances in the 2D-ES. If this were to happen, it might be wiser to take the 



approach of Cina and coworkers 
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- |91| to consider the evolution of the nuclear wavepackets 
for a few modes strongly coupled to the system, and maybe regard the rest of the modes as a 
)osonic bath. This possibility would require an exponential increase of experimental resources 



33|, so either partial or compressed sensing approaches 
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93| would be necessary. Alternatively. 



by going back to the time-domain picture provided by the authors in their previous work 



and applying novel concepts of QPT for initially correlated states 
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3^, 



a coarse grained and 



consistent tomographic protocol could be designed to address this problem. Finally, it might be 
worth considering additional nonlinear optical spectroscopic techniques such as considering the 



oscopic 



analysis of both rephasing and non-rephasing signals |59l. |97| . transient grating 2J|, pump probe 
or phase cycling of multipulse induced fluorescence |98| to investigate if they provide additional 
information for a more robust QPT. 

Albeit this article not being exhaustive, we hope to have convinced the reader that the QPT 
approach follows the spirit of MDOS in a very natural way. By systematically studying excited 
state dynamics as a quantum black box, an intriguing perspective on MDOS emerges that allows 
the use of tools designed in the QIP community in order to study excited state dynamics of 
condensed molecular systems. These possibilities will be the subject of future studies. 
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APPENDIX: Derivation of Eqs. ([6]), (|4]), and ([3]) 

Proof of Eq. (I6|).- Consider a system S interacting with a bath B. The total density matrix 
of the composite object is ptotai, whereas the reduced one for the system and the bath are p and 
Pb, respectively. Suppose that the total initial state is a tensor product of the form: 

Ptoia/(0) = p(0)®pb(0). (55) 
Where Pb(0) is assumed to be fixed at: 

PBiO) = J2Pp\ep){epl (56) 

with > 0, for every initial state p(0) of the system. 

At time T, the state of the composite object is simply a rotation of the initial state (it is a 
closed system): 

PtotaliT) = U{T)ptotal{0)U+{T). (57) 

Here, U{T) = Tie-'^o H,,tai{t')dt' ^ ^^le propagator for the entire object, where T is the time- 
ordering operator, and Htotai is given by: 

Htotai = Hs + Hb + HsB, (58) 

where Hs, Hb, Hsb are terms in the Hamiltonian that depend only on S, on B, or on degrees of 
freedom of both, respectively. Taking the trace of Eq. ( l57j) with respect to the states of B yields 
PiT): 

p{T) = Y,E^p{T)p{Q)Etf,{T) (59) 
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where: 



E^p{T) = ^^{e^\UiT)\e^), 



(60) 



is a Kraus operator and Eq. f l59p is often known as the operator sum representation 35|, |36|. By 
identifying: 



al3 



(61) 



we have shown the equivalence between Eq. (159|) and 



□ 



Proof of Eq. ((31). -As explained in Section 11 and in the previous proof, Eq. holds for 
any proper or improper state p(0) in the Liouville space of the system. With this in mind, Eq. 
([2]) is straightforward to prove. We want to enforce that Eq. preserves trace throughout the 
evolution in time T. If p(0) = \k){l\, then it must be that Tr{p(T)) = Tr{p{0)) = 5ki- Writing the 
elements of the initial density matrix as Pcd(O) = Sck^di, we immediately obtain the condition we 
want to prove: 



Tr{p{T)) = '^Sab'^XabcdSckSdh 



cd 



>kl 



Xaakl • 



(62) 
□ 



Proof of Eq. dH).- Again, we exploit the fact that Eq. (jH]) is true for any p(0). 
First, we consider a population as the initial state, i.e. p(0) = \k){k\ for any k. Note that p(0) 
is a proper density matrix, and hence Hermitian. Then, 



Pab{T) = Xabcd{T)5ck5dk 



cd 



Xabkk(T), 



(63) 
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PUT) = 2. \hacd 

dc 

= XbakkiT). 



(64) 



Since x(T) maps proper density matrices to proper density matrices, p{T) must also be Hermitian, 
i.e. pba{T) = p1^i^{T). It follows from Eqs. f l63l) and f l6^ that for any a, b, k: 

XbakkiT) = xlbkkiT). (65) 

This proves Eq. (jlj) for the case c = d. Next, we consider initial proper states of the form 
p(0) = ^{\k){k\ + \k){l\ + \l){k\ + 10(^1)- By repeating the exercise above, we get: 

^iXbakkiT) + Xbakl{T) + XbalkiT) + XballiT)) = ^{Xabkk{T) + XabkliT) + Xablk{T) + Xabll{T))* (66) 

Substituting the result Eq. (l65l) into Eq. (l66l) yields: 

Xbaki (T) + Xbalk (T) = xlbki (T) + xlbik (T) (67) 

This result almost proves Eq. (jlj) for the case c ^ d. In order to complete the proof, we analyze 
another class of initial proper states p(0) = l{\k}{k\+i\k){l\-i\l){k\ + \l){l\). The result from the 
same manipulations is that: 

iXbaki (T) - ixbalk (T) = -ixlbki (T) + ixluk (T) (68) 
Comparing Eqs. (1^ and (^7} yields: 



XbakliT) = xlbikiT) 

XbalkiT) = xlbkliT) (69) 



for any a, 6, /c, /. Eqs. (l65j) and fj69j) together yield Eq. 
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